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Abstract 

Inflation ends with the formation of a Bose condensate of inflatons. We show 
that in hybrid inflation models this condensate is typically unstable with respect 
to spatial perturbations and can fragment to condensate lumps. The case of D- 
term inflation is considered as an example and it is shown that fragmentation 
occurs if A > 0.2g, where A is the superpotential coupling and g is the U{V}fi 
gauge coupling. Condensate fragmentation can result in an effective enhance- 
ment of inflaton annihilations over decays as the main mode of reheating. In 
the case of D-term inflation models in which the Standard Model fields carry 
U{V)fi charges, if condensate fragmentation occurs then reheating is dominated 
by inflaton annihilations, typically resulting in the overproduction of thermal 
gravitinos. Fragmentation may also have important consequences for SUSY flat 
direction dynamics and for preheating. 
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1 Introduction 



A common feature of the cosmology of particle physics models is the formation of Bose 
condensates of scalar particles. Examples include axion condensates [||, condensates 
of squarks and sleptons along flat directions of the minimal supersymmetric (SUSY) 
standard model (MSSM) (Affleck-Dine condensates [Q, |3|, Q) and inflaton condensates 
which form at the end of inflation and whose decay is responsible for reheating the 
Universe 0. It is usually assumed that the scalar particles in the condensate are non- 
interacting, corresponding to coherent oscillations in a purely <fi 2 potential. However, 
in many cases this is not true. In the case of axions, deviation of the angular pseudo- 
Nambu Goldstone axion potential from a pure <ft 2 potential implies an attractive force 
between the axions which results in the growth of spatial perturbations and the for- 
mation of axion miniclusters ||. In the case of the Affleck-Dine condensate, deviation 
from the (p 2 potential, either due to the flattening of the potential above the messenger 
field mass (gauge- mediated SUSY breaking 0,3) or due to radiative corrections from 
gaugino loops (gravity- mediated SUSY breaking |], |10j) results the fragmentation of 
the condensate to form Q-balls || [], [II], O. Thus the conventional view of cosmo- 
logical condensates as being spatially homogeneous coherently oscillating scalar fields 
is not generally true. In particular, when the potential is 'natter-than-0 2 ', meaning 
Min(V '(0) / (p 2 ) is at 7^ (with V(0) = 0), the condensate is unstable with respect 
to spatial perturbations and fragments to non-topological solitons which we will refer 
to as condensate lumps. 

Here we consider the question of the stability of the inflaton condensate with respect 
to spatial perturbations and the consequences of its fragmentation. The most natural 



inflation models are hybrid inflation models [T2|], which, unlike the case of single-field 
inflation models, allow inflation to occur without requiring couplings to be very small. 
We will therefore focus on hybrid inflation models in the following. 

Although our results for inflaton condensate fragmentation can apply to hybrid 
inflation models in general, we will focus on the case of SUSY hybrid inflation models 



13], [LJ, pL5l , which will allow us to illustrate the general phenomenon of inflaton con- 
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densate fragmentation whilst applying the results to a case of considerable interest. 
SUSY hybrid inflation models are either of the F-term [HJ or D-term |T4| , |1(J type. 
The most interesting are the D-term models, which can evade the so-called 7/-problem 
i.e. the flatness of the inflaton potential in the presence of supergravity corrections 



HJ. We will therefore focus on D-term inflation, whilst presenting the results in a 
form that will allow them to be applied to other hybrid inflation models. 

Recently it has been shown that it is also possible for inflation to end via "tachyonic 
preheating" i.e. the rapid growth of spatial perturbations of the inflaton field in the 



presence of a tachyonic potential |17], [18|, |T||. The mode by which hybrid inflation 
ends (inflaton condensate fragmentation or tachyonic preheating) will be sensitive to 
the initial conditions at the phase transition ending hybrid inflation, in particular 
the rate of roll of the homogeneous scalar field relative to the rate of growth of the 
spatial perturbations. This requires a full analysis of the dynamics of the inflaton field, 



including the effect of radiative corrections to the inflaton potential [ED|. Since in this 
paper we wish to study the growth of spatial perturbations of a homogeneous hybrid 
inflation condensate in general, using D-term inflation as a particular example, we will 
assume throughout that a coherently oscillating scalar field condensate initially exists. 

The paper is organised as follows. In Section 2 we review the D-term hybrid 
inflation model. In Section 3 we discuss condensate instability in hybrid inflation 
models. In Section 4 we consider the evolution of spatial perturbations of a coherently 
oscillating condensate. In Section 5 we apply the results to the case of D-term inflation. 
In Section 6 we consider possible consequences of inflaton condensate fragmentation, 
in particular the enhancement of annihilations as a mode of reheating. In Section 7 we 
comment on the relationship between tachyonic preheating and inflaton condensate 
fragmentation. In Section 8 we present our conclusions. 

2 D-term Hybrid Inflation 



The superpotential of D-term inflation models is \14\ 



W = A5$ + $_ , (1) 
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resulting in a scalar potential 

2 

V = A 2 |S| 2 (|$+| 2 + |$_| 2 ) + A 2 |$+| 2 |$_| 2 + |- (|$+| 2 - |$_| 2 + e) 2 , (2) 

where Re(S) is the gauge singlet inflaton, Q± are fields with charges ±1 with respect 
to a Fayet-Illiopoulos U(l) gauge symmetry, U(1)fi, and £ > is the Fayet-Illiopoulos 
term. For \S\ > \S C \ = gy/l/\ the minimum of V(® + ,$-, \S\) is at $± = 0. With 
$± = 0, the tree-level S potential is flat with V = V = g 2 £ 2 /2 (^ 2 « 8.5 x 10 15 GeV 
from COBE normalization |15|). One- loop corrections result in a potential for S which 
causes S to slow-roll towards S = [13]. Once \S\ < \S C \, the minimum of the potential 



for a given value of \S\ is at $ + = and 



(In the following we may consider 5* and $_ to be real.) Thus the expectation value of 
$_ at the minimum of its potential is a function of the value of S. The mass squared 
terms along the S, $ + and directions as a function of \S\ and the $_ expectation 
value Eq. © are m 2 = A 2 |$_| 2 , m% + = A 2 |$_| 2 + 2A 2 |S| 2 and m 2 , = m\ = 2# 2 |$_| 2 , 
where m $ / is the mass at $_ ^ minimum Eq. (H) and A is the U(1)fi gauge boson. 



3 Condensate Instability in Hybrid Inflation Mod- 
els 

The dependence of the minimum of the $_ potential on the value of the S field is 
the reason for the instability of the inflaton condensate. Once \S\ < \S C \, S and $_ 
oscillate about the minimum of their potentials. Oscillations begin once rris > H . In 
the case of D-term inflation, this is satisfied once ISp/IS^I 2 = 1 — 47rg 2 £/3A 2 Mj, z = 
1 — 2 x 10~ 6 g 2 /A 2 (using mj = A|$_|, with |$_| as given by Eq. so S oscillations 
typically begin when \S\ is close to \S C \. The equation on motion for the inflaton, in 
terms of the conventionally normalized real scalar field s = ^/2Re(S), is 

V 2 

s + 3Hs- =rs = -A 2 s|$_| 2 , (4) 
a 2 
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where H = a/ a is the expansion rate and a is the scale factor. Suppose we consider the 
growth of a small spatial perturbation of s. As s decreases below s c , the mean value 
of the oscillating $_ field at a point in space will be approximately equal to the value 
at minimum of the $_ potential at that point in space, which depends of s(x, t). So 
if we average over the coherent oscillations of $_ about the minimum and replace $_ 
by the value $_ (s) at the minimum of its potential, the s equation of motion becomes 

V 2 A 4 s 3 
s + 3Hs -=rs^ -A 2 £s + — ^ . (5) 
a 2 2g 2 

Therefore the s scalar field and perturbations will evolve as if the s field had a potential 

\2t e 2 \4 e 4 

v.M « - • (6) 

This is a flatter than s 2 potential, corresponding to an attractive interaction amongst 



the s scalars and a negative pressure in the condensate [22|| . Therefore spatial per- 
turbations of the s condensate will grow, eventually becoming non-linear and resulting 
in fragmentation into condensate lumps || ^|, |1(J. The procedure of averaging over 
coherent oscillations of the field is well-defined if m^i is large compared with ms, 
which is true if is large compared with A, and we will focus on this case. In the 
case where one cannot first average over the $_ oscillations the combined dynamics 
of the S and $_ field will be more complicated. A particular case of this is F-term 
hybrid inflation, for which there is only a single coupling in the scalar potential such 
that the condition A = \plg is effectively satisfied [13, HJ. In this case there exists an 



exact solution of the scalar field equations such that the inflaton is described by an 
effective potential of the form as 2 — 6|s| 3 + as 4 (a, b, c > 0) ||24|| . 

Although we have derived V e ff(s) for the example of D-term inflation, we emphasize 
that a — s A attractive interaction is a generic feature of all hybrid inflation models for 
which we can average over the oscillations of the field terminating inflation prior to 
discussing the dynamics of the inflaton. Therefore our analysis may be readily applied 
to other hybrid inflation models. 
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4 Evolution of Perturbations 



We next consider the growth of spatial perturbations and the fragmentation of the 
inflaton condensate. The linear growth of perturbations has been discussed for a 
complex scalar field in the context of Q-ball formation in ||, using the approach of 
Here we adapt this approach to the case of a real scalar field in the expanding 
Universe. The equation of motion for a real scalar field $ is 

4 + !U rt_tf. = _«£W (7) 

a 2 <9$ 

We will assume throughout that V($) is a polynomial with V($) = V(— $). We define 
$ = (oo/«) 3//2 0, where a Q is the scale factor when the coherent oscillations begin. The 
equation of motion then becomes 

72 

a 2 



V 2 dUtf) 
-^<t>= JuT 



(8) 



where 

dU{<j)) ( a\ 3 / 2 dV($) 



+ A H , (9) 



&h = -1(h+Ih 2 ) . (10) 



where 

Arr = - 

2 V 2 

With = RSinQ, the equation of motion becomes 
(i?-i«i 2 -^-i?+ " v " t "" / )ffmf]+(-M]-2MH ~^%^"" +fl^^)Cosf] 



(ii) 

where Sj = d/dxi (i = 1,2,3). Multiplying the equation by SinQ and averaging over 
coherent oscillations gives 

A _aff_^ fl+ M = _«^a, (12) 

a 2 a 2 on 

where 
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and where < /(f2) >= ^- J 27r f(P) d£l denotes averaging over oscillations. (Ah is 

effectively constant on the timescale of coherent oscillations. In fact, since Ah is at 

most of the order of H 2 , the Ah term in general contributes a negligible correction to 

the mass squared term. Thus we take Ah = in the following.) Multiplying both 

sides by Cosfl and averaging gives 

t R- 2 djRdjV, _ 
R ~R~^ 2 a 2 



^ + 2-0--^-^ =-Q = 0. (14) 



In this we are assuming that R and fl do not vary much over the period of the oscil- 
lations. In practice we will be applying this method to the case of a — $ 4 interaction 
term in the potential. Therefore this method is accurate if the — $ 4 term is a small 
perturbation of the $ 2 term. (In the pure $ 2 limit R and Q are constant.) 

With R = R + 6R{x., t) and ft — Q(t) + <%7(x, t), the perturbation equations are 

V 2 / d 2 U \ 

5R - (l 2 5R - 2(l5(lR - =r5R = - - 8R (15) 

a 2 \ oR 2 1 

and 



R(t) 



R /?0 V 2 

^ + 2-^ + 2-^-2—^-^-^ = 0. (16) 
R R R 2 a 2 

Assuming the perturbations have the form 5R = 5R e s ^~ tlf " x , 8tt = 5fi e' s w~ tk " x 
[|, the perturbation equations become 

k 2 

(a + a 2 -Q 2 + — + U'' ff )5R = 2aR(l5tt (17) 
a 2 11 



and 



d + + — + 2— a Uf2 = 



a 



i? y \ R 2 R 



where a = S. Combining these gives a dispersion relation [25 , 
^ + « 2 -Q 2 + ^ + t/ e// J ^ + a 2 + - + 2- a j=2^^-^j . (19) 

In the case where there is no expansion the amplitude of oscillation is constant and so 
we have R = 0. A growing perturbation solution is then given by P5] , ^] 

k 2 ( Kff _ TT" \ 

a 2 = — X £ rsm 



and a = 0. In deriving this we have used (l 2 = U' e ^/R (from Eq. (|I2]) with constant 
R). This solution exists if k 2 /a 2 is less than k^^/a 2 = ^-|p- — ^C//^ • (I 11 obtaining 
Eq. ( |2T)| ) it is assumed that 16(k 2 / a 2 )U' e j f / R is small compared with (3U' e f j / R+U'^j j) 2 , 
which is satisfied for all k up to k max in the case where the — $ 4 potential term is small 
compared with the $ 2 term.) 

In the case with expansion we generally have to solve the equations of motion and 
perturbation equations numerically. However, for condensate fragmentation we will 
be mostly interested in the case where k 2 = k 2 ^, corresponding to the largest value 
of a at a given time and so the first perturbation mode to go non-linear. We will 
also be considering oscillation amplitudes such that the — $ 4 potential term is small 
compared with the $ 2 term. In this case R may be considered constant throughout. 
Then if at is non-zero, the solution Eq. ( |TT?D generalizes to 

k 2 ( U 'efJ Tj" \ 
a 2 I R U eff I 

« 2 = t ^-p (21) 

((47 - 1) if + K ff ) 

where a = a 2 + a and 7 = a 2 /a 2 . Since typically \a/a\ ~ H, we see that Eq. (f2~0|) 
will be approximately correct for the case k 2 = \^ max so long as a(k max ) > 
The solution of the perturbation equations is then 



H. 



5R « 8R exp (J adt^j e* kx (22) 

and 

5VI « 5n o exp ( f adt) e ik ' x . (23) 



We next apply the above to the case of the generic attractive — $ 4 interaction of 
hybrid inflation models, 



m 2 $ 2 ?7$ 4 

vm = - V ■ (24) 



From Eq. (|T§, 



«/„,(*> _ 2 / v _ #4 = 2 /™w _ 3^ 



<9i? i?\ V a / / R \ 2 8 V a 
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where we have used < Sin fl >= 1/2 and < Sin Hi >= 3/8. Thus 

u ^ = — 2 ifW • (26) 

Thus from Eq. ( |2Ti| ) we find 

a \ 3 / 2 f3r]R 2 \ l/2 / 9//.R 2 /a \ 3 \~ 1/2 |k| 



«= - SH i-Sr - -■ (27) 



ay \ 8m 2 / \ 8m 2 



Eq. (|27|) is strictly valid only if the — $ 4 term is a small perturbation of the $ 2 
term. In this case R will be essentially constant (equal to its initial value R Q ) and 

~~ (if) ) a PP rox i ma tely equal to 1, conditions which we will assume 

to be satisfied in the following. 

The growth of the perturbations is then given by, 

3r]R 2 \ 1/2 |k| / 2 \ /a \ 3 / 2 -™ //a\ 5 /2-" 



6R 6R 
— « — exp 
it rt 



m 2 J aH a V 5 — 2nJ \ a J \\a 



(28) 



In this we have used H = H (a /a) n , where n will be between and 3/2 as the <3> 
oscillations develop from the end of inflation to an approximately $ 2 potential. The 
condition for fragmentation to occur is that 5R/R > 1. The largest growth at a given 
time corresponds to the mode k max , where 



a 2 ~ \ R " eff ~ 2' 



a 



-Ue ff =~VR 2 ~ ■ (29) 



This determines the radius of the condensate lumps when the condensate fragments, 

n, 

vra /2\V2/ a x3/2 n 

ri ~K^\~\3j UJ ^rW m 

This is really the initial radius of the lumps immediately after fragmentation, and 
the lump will subsequently relax to its stable configuration, in which the attractive 
potential term is balanced by the gradient term in the equation of motion. However, 
the radius of the stable configuration is similar to that of the initial lump, since 
for a stable configuration of the form <f)(r,t) « <ft(r)Sin(mt), the equation for <f)(r) is 
d 2 cj)(r)/dr 2 +(2/r)d(j)(r)dr « Sin^im^dSV/dcf), where V = m 2 (j) 2 /2+5V . For a stable 
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lump of radius r s and field amplitude we therefore expect that r 2 « <p\d5V / d(f)\~ l 
(the left hand side of the <f>(r) stable lump equation being ~ <p/r 2 ). For 5V = —rjtfi/A 
and (j) ~ R this implies that r s m (rjR 2 )^ 1 ^ 2 . 

To give a condition for inflaton condensate fragmentation, we use k = k maa; . The 
condition for fragmentation to occur is then 

2 

As a increases, the left hand side of Eq. fl3T|) is maximized for ^ = ( 2 (3- n ) ) 5 ~ 2ra anc ^ 
then decreases. Condensate fragmentation in hybrid inflation models must therefore 
occur soon after coherent oscillations begin if it is to occur at all. 

Assuming that the — <3> 4 contribution to the potential is small, the coherent oscil- 
lations are approximately $ 2 and so n = 3/2, such that maximum growth occurs at 
— = ~. The condition for condensate fragmentation is then 

a 3 ° 



1 r,Ri 



> 



2^3 mH, 

where H D is calculated using for the energy density p = m 2 R 2 /2. 



(3 , (32) 



5 Application to D-term Inflation 

For the case of D-term inflation, m 2 = A 2 £ and rj = X 4 /2g 2 . In order that the 
— $ 4 contribution to the potential is small, we choose R Q = s c /2, such that R 2 = 
g 2 £/2A 2 . (The perturbations will start to grow as soon as coherent oscillations begin 
with amplitude s ~ s c . However, as we cannot apply our method to calculate the 
growth for oscillation amplitudes between s c and s c /2, we are in fact underestimating 
the total growth.) Thus from Eq. (|32|) the condition for condensate fragmentation 
becomes 

^ > „ 0.014 . (33) 

9 P Mpi 

In order to complete the fragmentation condition we need the value of (3. The seed 
perturbations of the inflaton field are expected to come from quantum de Sitter 
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fluctuations of the scalar field during inflation. Modes with wavenumber k large 
compared w ith H will be excited by the increase of the inflaton mass at the end 
of inflation from approximately zero to ms ~ A^ 1 / 2 . The largest wavenumber ex- 



cited typically corresponds to k m pa m s , with amplitude 8s ~ m s /(27r) [^3|. Since 



|k ma:r |/a ~ y 3/2A£ 1//2 ~ k m at a pa a , seed perturbations of wavenumber k max /a will 
exist after inflation ends, with 

where we have used 5R /R ~ 5s/ s c . So with A, g > 0.1, /3 ~ 5 — 10 is a typical value. 
Therefore, using (3 = 10, we find that the condition for fragmentation to occur is 

A > 0.2^ . (35) 

Since this neglects the growth of perturbations for oscillation amplitudes between s c 
and s c /2, the true lower bound on A for fragmentation is likely to be smaller. The 
radius of the condensate lump relative to the horizon radius when the condensate 
fragments is then 

7r k (-J fkr 37x10 !■ (36) 

The condition for the approximations leading to Eq. ( p2|) to be consistent, a(k max ) > 
^ = 5H/2 at ^ = | (where we assume H oc a -3 / 2 for s < s c /2 and we have used 
a oc a~ 5//2 from Eq. fl27D ), is satisfied if 



A >4 0^, QQ7 (37) 
g M Pl 

So if the fragmentation condition Eq. ( |35D is satisfied then the approximations are 
consistent. 

Therefore in D-term inflation models condensate fragmentation is likely to occur 
if A ^ 0.2g. Since the true lower bound on A for fragmentation is likely to be smaller, 
it is probable that condensate fragmentation will occur if A and g both take values in 
the natural range 0.1-1. 
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6 Enhancement of Inflaton Annihilations By Con- 
densate Fragmentation 

We next consider the possible consequences of inflaton condenate fragmentation. One 
potentially important consequence is that inflaton annihilations are effectively en- 
hanced compared with the case of a homogeneous condensate and may dominate over 
decays as the primary mode of reheating. 

We first show that once the condensate fragments the particles within the conden- 
sate lump decouple from the expansion of the Universe, such that the number density 
and field amplitude inside the condensate lumps remains constant. This will be true 
if the force on the particles due to scalar interactions is greater than the gravitational 
force responsible for slowing the expansion of the Universe. Suppose the energy den- 
sity is dominated by a pressureless homogeneous energy density p. (For the case of 
inflaton condensate fragmentation this is not strictly true, since Sp/p ~ 1 when the 
condensate fragments. However, it is useful to make this assumption in order to ob- 
tain an expression that can be applied to scalar field models in general.) Suppose then 
we consider a spherical lump of radius r. The gravitational acceleration acting on a 
particle at the surface of the lump is then 



If the force due to the attractive scalar interaction produces a smaller acceleration 
than this, the particles will follow the expansion of the Universe, otherwise they will 
decouple from expansion. To estimate the force due to the scalar interaction consider 
a sphere of radius r and with a fixed number of scalar particles N, 



where the number density of scalars is n = m<j) 2 /2 and where for simplicity we have 
considered a constant amplitude for the coherently oscillating field inside the sphere. 
This gives as a function of r. With V{4>) = ' nL §- + SV, the total potential energy of 




(39) 



If 



the sphere is 

Thus the force due to the scalars is 

F . = _!?5- = _f^-^. (4!) 
dr V mr y 

For the case 5V = — ?70 4 /4, the force is F s = —irrj^r 2 . The condition that the 
acceleration due to scalar attraction is larger than the gravitational acceleration is 
then 

(42) 

Since condensate fragmentation occurs soon after coherent oscillations begin (at dV/d(j) • 
0), we have 2 ~ m 2 /r]. With r rs r; (with a a Q in Eq. (|30|)) Eq. becomes 

to 3 /S\ x / 2 /7 2 

>^ Sr- (43) 



r] 3 / 2 R \2J 2tt 

Note that if this condition is satisfied then from Eq. (^) we have \r/r\ ^ H 2 , which can 
be rewritten as \f\H~~ 1 > Hr = r. Therefore the attractive force between the scalars 
will bring the expansion of the lump to a halt within an expansion time, 5t ~ H^ 1 . 
For the case of D-term inflation, Eq. (|I3| ) is satisfied if 

A<(^) 1/2 «4xlOS (44) 

which is generally strongly satisfied. Thus the scalars in the condensate lumps decouple 
from the gravitational expansion. 

To see how condensate fragmentation effectively enhances annihilations, consider 
an interaction A^SplQiJ 2 between inflatons and light scalars Qi (in D-term inflation Qi 
will correspond to MSSM scalars 0]). We will calculate the annihilation rate assuming 
no Bose enhancement (i.e. no parametric resonant decay of the condensate), in which 
case the annihilation rate simply corresponds to the perturbative annihilation rate 
of the scalars in the condensate. The average of the annihilation cross-section times 
relative particle velocity in the v — > limit is JM] 

A 4 

< OV >ann= ~. K ■ (45) 
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With the s number density in the condensate given by n = mss 2 /2, the annihilation 
rate of scalars in the condensate is then 

A 4 s 2 

Tann = U < <JV > a nn= T7^" • (46) 

12S7rm s 

The condition for scalars to annihilate is then T ann > H = H (a /a) 3 ^ 2 . We can now 
see why annihilation is effectively enhanced when the condensate fragments. In the 
case of a homogeneous condensate, the scalars are freely expanding and s 2 oc a~ 3 , so 
the annihilation rate drops more rapidly than the expansion rate as the scale factor a 
increases. Thus unless annihilations are effective immediately after the end of inflation, 
they will never be significant and inflaton decays will be the main mode of reheating. 
However, if the inflaton condensate fragments soon after oscillations begin then the 
value of the amplitude s inside the lumps is constant and so T ann /H oc a 3 / 2 increases 
as a increases. Therefore as the Universe expands annihilations will eventually occur 
and may be the dominant mode of reheating. The reheating temperature due to 
annihilations is then 



where we have used s ~ s c inside the lumps and where g(Td) is the number of degrees 
of freedom in thermal equilibrium ||. 

We next apply this to the specific case of D-term inflation. If the MSSM fields also 
carry U(1)fi charges, there is an interaction with the scalars of the MSSM of the form 
A 2 |<3j| 2 |5l 2 coming from integrating out the U(1)fi gauge fields W&, where A is the 



superpotential coupling from Eq. (0). The reheating temperature is then 
Mir) \4n 3 g(T d )) Pl \g(T R 



gX^^M 1 ^ w 5.5 x 10 15 ( ) gX 1 ' 2 GeV . 

(4? 



The upper bound from requirng that gravitinos are not generated excessively by ther- 
mal scattering is Tr < 10 8-9 GeV p7|. Thus we require 



2A 1/2 < 1.8 xlO- 7 f4^V /4 (77^7) • (49) 



100 / V10 9 GeV 
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The smallest value of A for which fragmentation occurs is A ~ 0.2g. In this case if 
condensate fragmentation occurs then from Eq. ( f49|) the U(1)fi gauge coupling must 
satisfy g < 7 x 10~ 5 . (This upper bound will be even stronger in the presence of 
parametric resonance.) This will not be satisfied if, for example, the U{1)fi gauge 
coupling has the typical magnitude g ~ 1 of the MSSM gauge couplings. Thus con- 
densate fragmentation is typically not compatible with D-term inflation if the inflatons 
can annihilate to MSSM fields. In order to have D-term inflation consistent with the 
absence of thermal gravitinos we must either eliminate inflaton annihilations, which 
requires that the MSSM fields do not carry U(1)fi charges, or eliminate inflaton con- 
densate fragmentation, which requires that A < 0.2g. 

7 Inflaton Condensate Fragmentation and Tachy- 
onic Preheating 

Throughout the preceeding discussion we have assumed that a homogeneous inflaton 
condensate forms at the end of inflation. However, it has been shown that spatial 
perturbations of the inflaton field may grow and become non-linear much more rapidly 
(before homogeneous oscillations are established) in a process known as tachyonic 
preheating [17], In this case the final state is composed of colliding scalar field 
waves |T7|, |T^]. Recently it has also been shown |H5| that condensate lumps occur 
in tachyonic preheating (called "oscillating hot spots" in [|l9j)- Tachyonic preheating 
typically occurs in less than the time for a single coherent oscillation [[ij]], making it 
impossible to average over a coherently oscillating inflaton field. Inflaton condensate 
fragmentation and tachyonic preheating may be regarded as different manifestations 
of a general phenomena, namely the instability of the inflaton field with respect to 
spatial perturbations in hybrid inflation models. 

The question of whether inflaton condensate fragmentation or tachyonic preheating 
occurs at the end of hybrid inflation will depend upon the state of the field at the time 
when the inflaton reaches s c , in particular the rate at which the homogeneous field 
s(t) is rolling (where s — s(t) + 5s(x, t)) relative to the rate of tachyonic growth of the 
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spatial perturbations 5s(x, t) at s < s c . If the homogeneous field can catch up with 
the spatial perturbations (which cross s c before the homogeneous field) before there is 
significant growth of the perturbations due to the tachyonic mass term at s < s c , then 
there will be no tachyonic preheating. This depends crucially on the rate of rolling of 
the homogeneous field at s c , which in turn requires that the full inflaton potential with 
radiative corrections be considered. For the case of D-term inflation a full analysis has 
yet to be done EC . 



8 Conclusions 

We have shown that it is possible for the inflaton condensate in hybrid inflation models 
to fragment to condensate lumps. The inflaton condensate is in general unstable, but 
the instability reduces as the Universe expands, requiring that fragmentation occurs 
shortly after the end of inflation. The state of the Universe after fragmentation, with 
the energy density concentrated inside inflaton condensate lumps, is quite different 
from the conventional post- inflation scenario of a homogeneous inflaton condensate. 
One consequence of inflaton condensate fragmentation is that inflaton annihilations 
will be effectively enhanced relative to the case of a homogeneous inflaton condensate. 
In the case of D-term inflation models, which we have used as a specific example of 
hybrid inflation models in our discussion, if condensate formation and fragmentation 
occurs then the enhancement of inflaton annihilations implies that the reheating tem- 
perature is typically large compared with the thermal gravitino upper bound. Thus 
either inflaton annihilations must be suppressed, which requires that the MSSM fields 
do not carry U(1)fi gauge charges, or fragmentation must not occur, which requires 
that A < 0.2#. 

If inflaton annihilations do not occur, or if we consider a non-SUSY model for 
which there is no thermal gravitino bound on the reheating temperature, then inflaton 
condensate fragmentation can safely occur. This may have interesting consequences for 
the inflaton dominated period following inflation. For example, the fact that the energy 
density of the Universe is now concentrated in condensate lumps could greatly alter the 
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dynamics of SUSY flat direction scalar fields, which acquire a mass from the SUSY 
breaking inflaton energy density in the case of a homogeneous inflaton condensate 
0]. Another potentially important effect would be for parametric resonant decay of 
the inflaton and preheating. In the case of a homogeneous condensate, parametric 
resonance turns off as the Universe expands and the inflaton oscillation amplitude 
decreases |28[ |29| . However, if the condensate fragments then, just as in our discussion 
of perturbative annihilations, the oscillation amplitude will be frozen inside the lumps 
and so parametric resonant decay should continue without stopping, resulting in more 
efficient preheating. 

Growth of spatial perturbations of the inflaton field in hybrid inflation models has 
also been demonstrated in the context of tachyonic preheating. In general, the mode 
by which hybrid inflation ends is likely to be sensitive to both the model paramaters 
and the initial conditions at the end of inflation. A detailed numerical investigation 
will therefore be necessary in order to clarify how hybrid inflation ends in a given 



model | 20| . However, it should be emphasized that condensate lumps are a feature 



of both inflaton condensate fragmentation and tachyonic preheating fL9| . Therefore 



our discussion of the cosmology of condensate lumps should apply to the tachyonic 
preheating case also. 

It is important to emphasize that inflaton condensate fragmentation is a natural 
possibility in all hybrid inflation models and that it may have consequences for cos- 
mology beyond what has been discussed here. A detailed understanding of inflaton 
condensate fragmentation will therefore be necessary in order to fully understand the 
cosmology of hybrid inflation models. 
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